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Abstract. We present the exact solutions of various directed walk models of polymers 
confined to a slit and interacting with the walls of the slit via an attractive potential. 
We consider three geometric constraints on the ends of the polymer and concentrate 
on the long chain limit. Apart from the general interest in the effect of geometrical 
confinement this can be viewed as a two-dimensional model of steric stabilization and 
sensitized flocculation of colloidal dispersions. We demonstrate that the large width 
limit admits a phase diagram that is markedly different from the one found in a half- 
plane geometry, even when the polymer is constrained to be fixed at both ends on 
one wall. We are not able to find a closed form solution for the free energy for finite 
width, at all values of the interaction parameters, but we can calculate the asymptotic 
behaviour for large widths everywhere in the phase plane. This allows us to find the 
force between the walls induced by the polymer and hence the regions of the plane 
where either steric stabilization or sensitized flocculation would occur. 
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Short title: A directed walk model of a long chain polymer in a slit with attractive 
walls 
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1. Introduction 

When a long linear polymer molecule in dilute solution is confined between two parallel 
plates the polymer loses configurational entropy and exerts a repulsive force on the 
confining plates. If the monomers are attracted to one of the two confining plates, 
the force is still repulsive, but less is magnitude because the polymer can adsorb on 
one plate and therefore does not extend as far into solution. The effect of the second 
confining plate is then less and the loss of configurational entropy is less. If the monomers 
are attracted to both plates the net force can be attractive at large distances and 
repulsive at smaller distances. At large distances the energy term is dominant while 
at smaller distances the entropy loss is dominant. These phenomena are related to the 
stabilization of colloidal dispersions by adsorbed polymers (steric stabilization) and the 
destabilization when the polymer can adsorb on surfaces of different colloidal particles 
(sensitized flocculation). 

One would hope to be able to say something about this problem for a relatively 
realistic model of a polymer in a good solvent such as a self-avoiding walk and the 
self-avoiding walk case has been considered by several authors (see for instance Wall et 
al 1977 and 1978, Hammersley and Whittington 1985). When the self- avoiding walk is 
simply confined between two parallel lines or planes, and does not otherwise interact 
with the confining lines or planes, a number of results are available. Let c n (w) be the 
number of self-avoiding walks on the simple cubic lattice Z 3 , starting at the origin and 
confined to have the z-coordinate of each vertex in the slab < z < w. Then it is 
known that the limit 

lim n -1 log c n (w) = k(w) (1.1) 

exists. If c n is the number of n-edge self-avoiding walks with no geometric constraint 
then the connective constant of the lattice is given by 

lim n _1 logc„ = k. (1.2) 

n— +oo 

It is known that 

(i) k(w) is monotone increasing in w, 

(ii) lim^^oo k(w) = k, and 

(iii) k(w) is a concave function of w. 

In two dimensions, the exact values of k(w) are known for small w (Wall et al 1977). 
When the walk interacts with the confining lines or planes almost nothing is known 
rigorously but the problem has been studied by exact enumeration methods (Middlemiss 
et al 1976). 

To make progress with the situation in which the walk interacts with the confining 
lines or planes one has to turn to simpler models. In a classic paper DiMarzio and 
Rubin (1971) considered a random walk model on a regular lattice where the random 
walk is confined between two parallel lines or planes, and interacts with one or both of 
the confining surfaces. They see both attractive and repulsive regimes in their model, 
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though they confine themselves to the case where either the polymer interacts with only 
one surface or equally with both surfaces. 

In this paper we study a directed version of the self-avoiding walk model. We 
consider directed self-avoiding walks on the square lattice, confined between two lines 
(y — and y = w). We consider the cases where the walk starts and ends in y — 
(confined Dyck paths, or loops), starts in y — and ends in y = w (bridges) and starts 
in y = but has no condition on the other endpoint beyond the geometrical constraint 
(tails). Our model is related to that of DiMarzio and Rubin (1971) but we consider the 
more general situation where the interaction with the two confining lines or planes can 
be different. The techniques which we use are quite different from those of DiMarzio 
and Rubin. 

We derive results for the behaviour of these three classes of directed walks as a 
function of the interaction of the vertices with the two confining lines, and the width 
w. In particular, we find the generating functions exactly for each of the three models. 
We find the free energy and hence the phase diagram in the limit of large slit width for 
polymers much larger than the slit width. We demonstrate that this phase diagram is 
not the one obtained by considering the limit of large slit width for any finite polymer. 
We calculate the asymptotics for the free energy in wide slits and so deduce the forces 
between the walls for large slit widths. This allows us to give a "force diagram" showing 
the regions in which the attractive, repulsive, short-ranged and long-ranged forces act. 

2. The model 

The directed walk model that we consider is closely related to Dyck paths (see Stanton 
and White 1986 or Deutsch 1999). These classical objects are closely related to Ballot 
paths which were studied as a voting problem (Andre 1887 and Bertrand 1887). Dyck 
paths are directed walks on Z 2 starting at (0, 0) and ending on the line y = 0, which 
have no vertices with negative ^-coordinates, and which have steps in the (1, 1) and 
(1, —1) directions. We impose the additional geometrical constraint that the paths lie 
in the slit of width w defined by the lines y = and y = w. We refer to Dyck paths 
that satisfy this slit constraint as loops (see figure HJ. We also consider two similar sets 
of paths; bridges (see figure EJ) and tails (see figure EJ). A bridge is a directed path lying 
in the slit which has its zeroth vertex at (0, 0) and last vertex in y = w. A tail is a 
directed path in the slit whose zeroth vertex is at (0, 0) and whose last vertex may lie 
at any < y < w. 

Let C w , B w and T w be the sets of loops, bridges and tails, respectively, in the slit of 
width w. We define the generating functions of these paths as follows: 




(2.1) 




(2.2) 
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Figure 1. An example of a directed path which is a loop: both ends of the walk are 
fixed to be on the bottom wall. 




Figure 2. An example of a directed path which is a bridge: one end of the walk is 
fixed to be on the bottom wall (y = 0) while the other end is fixed to be on the top wall 

(y = w). 




Figure 3. An example of a directed path which is a tail: one end of the walk is fixed 
to be on the bottom wall (y = 0) while the other end is not constrained. 



T w (z,a,b) = J2z n(p) a u(p) b v(p) , (2.3) 
peT w 

where n(p),u(p) and v (p) are the number of edges in the path p, the number of vertices 
in the line y = (excluding the zeroth vertex) and the number of vertices in the line 
y = w, respectively. 

Let and be the sets of loops, bridges and tails of fixed length n in the 

slit of width w. The partition function of loops is defined as 

Z l n oop (w; a, b) = a u(p) b v{p) (2.4) 
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with the partition functions for bridges and tails defined analogously. Hence the 
generating functions are related to the partition functions in the standard way, eg for 
loops we have 

L w (z,a,b) = J2^Z l n ° op (w,a,b). (2.5) 

n 

We define the reduced free energy K loop (w; a, b) for loops for fixed finite w as 

K loop (w; a, b) = lim n" 1 log Z l n oop {w; a, b) (2.6) 

n— »oo 

with the reduced free energies for bridges K brid9e (w; a,b) and tails K tml {w;a,b) defined 
analogously. 

Consider the singularities z l ° op (w; a,b), z brid9e (w; a,b) and zl ml (w ; a, b) of the 
generating functions L w (z,a,b), B w (z,a,b) and T w (z,a,b), respectively, closest to the 
origin and on the positive real axis, known as the critical points. Given that the radii 
of convergence of the generating functions are finite, which we shall demonstrate, and 
since the partition functions are positive, the critical points exist and are equal in value 
to the radii of convergence. Hence the free energies exist and one can relate the critical 
points to the free energies for each type of configuration, (that is loops, bridges and 
tails) as 

n type (w; a,b) = - log z t9pe (w; a, b). (2.7) 

Theorem 2.1. Loops, bridges and tails have the same limiting free energy at every fixed 
finite w. 

Proof. We sketch the proof: Fix w at some finite positive integer, and let Z 1 "^, ^bridge 
and Z^ d be the partition functions of loops, bridges and tails with n edges in the strip 
of width w. By appending w edges to loops and bridges we have the inequalities 

bZ l r < Z^J e and aZ b n " d ^ < Z l n °Z- (2.8) 

Hence abZ\°° v < aZ h ^+^ e < Z l °^ w . Taking logarithms, dividing by n and taking the 
limit n — > oo shows that the limiting free energies are equal (since the limits exist). 

Since every loop is a tail Z l °° v < Z^ 1 , and by appending w (if n and w have the 
same parity) orw-1 edges (if n and w have different parity) to tails we have 

aZ tau < z io Z or aZ^ l <Z l n Z^ (2.9) 

depending on the parity of n + w. By a similar sandwiching argument the limiting free 
energies of loops and tails, and hence bridges, are all equal. □ 
Since the free energies are the same for the three models we define k(w; a, b) as 

k(w; a, b) = K loop {w; a, b) = K brid9e {w; a, b) = K tm \w; a, b) (2.10) 

with z c (w; a, b) being 

z c {w; a, b) = z l c oop {w; a, b) = z b " d9e {w; a, b) = z l f[w] a, b) (2.11) 

so that 

k(w; a, b) = — \ogz c (w; a, b). (2-12) 
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3. The half-plane 

For loops and tails we first consider the case in which w — > oo, which reduces the 
problem to the adsorption of paths to a wall in the half-plane. We begin by noting that 
once w > n for any finite length walk there can no longer be any visits to the top surface 
so 

C n w+1 = Cl and T: +1 = T: for all w > n. (3.1) 

Hence we define the sets C^ p and T£ as 

^ P = C+i and T-=T: +l . (3.2) 

The limit w — > oo can be therefore be taken explicitly. Also, as a consequence of the 
above, for all w > n we have v(p) = for any p e £™ and also for any p e T™. Hence 
the partition function of loops in the half-plane can be defined as 

Z l n oop ' hp {a) = ( 3 - 3 ) 

with the partition function of tails in the half-plane defined analogously. We define the 
generating function of loops in the half plane via the partition function as 

L(z,a)=J2 znZ n 0P ' hP ( a ) ( 3 - 4 ) 

n 

with the generating function of tails T(z, a) in the half-plane defined analogously. Note 
that the above limit does not exist for bridges since for any fixed walk of length n we 
have w < n. In an analogous way to the slit we define the reduced free energy in the 
half-plane for loops (and tails analogously) as 

i) = lim n" 1 \ogZ l n oop ' hp (w; a). (3.5) 



loop, hp I 



K 

We note that in defining these free energies for the half-plane the thermodynamic limit 
n — > oo is taken after the limit w — > oo: we shall return to this order of limits later. 

While the half-plane solution is well known (Brak et al 1998, 2001 or Janse van 
Rensburg 2000) it is useful to summarise the main results for comparison with our new 
results for slits. One may factor both loops and tails in the half plane by considering the 
second last point at which a loop touches the line y = or the last point at which a tail 
touches the surface. This leads to the functional equations L(z, a) = l+az 2 L(z, a)L(z, 1) 
and T{z, a) = L(z, a) + zL{z, a)T(z, 1) for the generating functions for loops and tails. 
These may be easily solved to give 

2 

Liz, a) — . and (3.6) 

V ; 2 - a + ay/ 1 - 4z 2 V ; 

T(z, a) = A Z (3.7) 

y J (a- 2 -ay/1- 4^ 2 )(1 -2z-y/l- ^z 2 ) V ' 

The generating functions are singular when the square root term is zero and when the 
denominator is zero. These conditions give the locations of the critical points, z l ° op ' hp (a) 
and z* mZ ' fep (a) and so also give the free energies. The values are equal for the two models 
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and so we define z^ p (a) = z l ° op ' hp (a) = z t ™ l,hp {a), and similarly for the free energy we 
have 



K hp {a) = K loop ' hp {a) = K tail ' hp (a) = - log z^{a) (3.8) 



with 



f log(2) a<2 

The second derivative of K hp (a) is discontinuous at a = 2 and so there is a second-order 
phase transition at this point. The density of visits in the thermodynamic limit, defined 

as 

^ (o) = „,„ M = *£M (3.10) 
n->oo n d log a 

is an order parameter for the transition and has corresponding exponent (3 = 1. 

4. Exact solution for the generating functions in finite width strips 

The solution of the strip problem proceeds in an altogether different manner to the 
half-plane situation. This begins by using an argument that builds up configurations 
(uniquely) in a strip of width w+ 1 from configurations in a strip of width to. In this way 
a recurrence-functional equation is constructed rather than a simple functional equation. 
Consider configurations of any type — loops, bridges or tails — in a strip of width w 
(see figure EJ), and focus on the vertices touching the top wall: call these top vertices. 
These vertices contribute a factor b to the Boltzmann weight of the configuration. Now 
consider a zig-zag path (see figure EJ), which is defined as a path of any even length, or 
one of length zero, in a strip of width 1. The generating function of zig-zag paths is 
1/(1 — bz 2 ). Replace each of the top vertices in the configuration by any zig-zag path. 
Since one could choose a single vertex as the zig-zag path all configurations that fit in a 
strip of width w are reproduced. Also, the addition of any non-zero length path at any 
top vertex will result in a new configuration of width w + 1 and no more. The inverse 
process is also well defined and so we can write recurrence-functional equations for each 
of the generating functions. The generating function L w (z,a,b) satisfies the following 
functional recurrence: 

Li(z,a,b) =- — -— ; (4.1) 
1 — abz z 

L w+1 (z, a, b) = L w ^z, a, - • ( 4 - 2 ) 

The generating function of bridges, B w (z, a, b), satisfies a similar functional recurrence: 

bz 



B w+ i(z, a, b) = bz B w [z,a, — ^ ) . (4.4) 
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Figure 4. The construction of configurations in strip width w + 1 from configurations 
in strip width w is illustrated. Every vertex touching the upper wall in the strip of 
width w can be replaced by a zig-zag path as shown. 



We note that the zeroth vertex of the path is weighted 1 and that L w (z,a,b) counts 
the walk consisting of a single vertex. To compute the generating function of tails we 
introduce a new variable c such that a tail whose last vertex is in the line y = w is 
weighted by c. Using this variable we arrive at the following functional recurrence: 

f^a^c) ={^J; (4-5) 

f u (z 1 o 1 6,c)=f u Uo,^,^j. (4.6) 

The relevant tail generating function is then T w (z, a, b) = T w (z, a, b, b). 

Using these results it is easy to prove by induction that the generating functions 
for any finite w must be ratios of polynomials in z. In fact using induction it is possible 
to prove that the generating functions L w (z, a, b), B w (z, a, b) and T w (z, a, b, c) have the 
following forms: 

L ^ h) =pjz^by (47) 

bz w 

B w (z,a,b) =— -, (4.8) 

P w (z,a,b) 

T w (z,aAc) = Q jl Z ' bj % (4.9) 
P w {z,a,b) 

where the P w (z, a, b) and Q w (z, b, c) are polynomials. By simply iterating the recurrence- 
functional equations and using, for example the GFUN package of MAPLE™ or by 
combinatorial means (Flajolet 1980 and Viennot 1985), one can find the recurrences 
Pw+i = Pw~ z 2 P w -i and Q w+ i = (1 + z)Q w - z(l + z)Q w _ 1 + z 3 Q w _ 2 . The solution of 
these, using yet another generating variable t over width w, gives 

G p (z, a, b, t ) = V P w (z, a, b)t w = t x ^ J \ >— (4.10) 

L — ' 1 — t + z z t z 

w>l 
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and 

/ , \ ^ / , x ,„ (1 + cz ) — zt(c + bz) + cz 3 t 2 /,-,-, x 

G g (z, a, b, c, t) = Qw{z, b, c)t w = t x 1 (i _ ^(i _ t + ^2) • ( 4 - U ) 

Alternately, one can find from substitution into the functional-recurrences f!4.2|) and 
flUJ) that 

P w+1 {z, a, b) = (1 - bz 2 )P w {z, a, Y^i) ( 4 - 12 ) 

and 

1 1 + cz 

Q w+1 (z, b, c) = (1 - bz 2 )Q w (z, — ^ , j-^j). (4.13) 

By multiplying both sides by £ w and summing over w, one obtains functional equations 
for G p and G g , which presumably can be solved directly. We note that P w and Q w 
are polynomials related to Chebyshev polynomials of the second kind (see for instance 
Abramowitz and Stegun 1972 or Szego 1975) as can be seen by comparing generating 
functions. Other techniques such as transfer matrices (Brak et al 1999), constant 
term (Brak et al 1998) and heaps (Bousquet-Melou and Rechnitzer 2002) would also 
undoubtedly work on these problems. 



5. Location of the generating functions singularities 

5.1. Transformation of the generating variable 

We recall that the limiting free energy of a model is determined by the dominant 
singularities of its generating function. Since the generating functions of these three 
models are all rational with the same denominator, P w (z,a,b), it is the smallest real 
positive zero of this polynomial that determines the limiting free energy. Theorem 12.11 
shows that this zero (since the only zero of the numerator for bridges is itself) is not 
cancelled by a zero of the numerator. Below we find the limiting free energy by studying 
the zeros of P w (z, a, b). Given w we would like to find an expression for z c as a function 
of a, b. 

Since the generating function G p (t) for the P w polynomials is a rational function of 
t we can write it in partial fraction form with respect to t, and then expand each simple 
rational piece to find P w as a function of w. The decomposition of the right-hand side 
of ()4.10|) involves expressions containing surds which becomes quite messy. It is easier 
first to make a substitution that greatly simplifies the partial fraction decomposition 
and the subsequent analysis. This is chosen so as to simplify the denominator of G p (t). 
If we set 

1+q y ' 
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the generating function of the polynomials then simplifies to 

[1 + q - aq)(l + q - bq) 



G P {^/q/(l + q),a,b,t) = ab-a-b-, _ q \( 1 + q _ t \ 

'l + q-a)(l + q-b) 



(l-q)(l+q-tq) 

where there are simply linear factors in the denominator. The partial fraction expansion 
has been taken in 1)5. 2|) . In our discussion of the singularities of P w we will start with 
this g-form. We can convert our results back later by making the inverse substitution 

1 - 2z 2 - VT^4^ 
9= ^ • ( 5 - 3 ) 

The above partial fraction form implies that 

P(/a/(l + a)ab)- (1 + 9 ~ + 9 ~ &g) - (1 + Q - q)(1 + Q - % ^ f5 4) 
and when a = this simplifies to 

Pf /a/(l + a) 6)- + + (55) 

We also have 

q h+™ (1 -6 + g) 



c 



(l-g) (l + qf 
1 + q — bq 



+ 



gf (1 - c - (b - 2 c) y/q + (1 - c) q) 



(5.6) 

(i-v^) (i+?r 

We can therefore find the generating functions for loops, bridges and tails. The 
generating function L w (y/q/ (1 + q), a, b) can be written as 

T _ (l + q)[(l + q-bq)-(l + q-b)q™} 



"' (1 + q ~ aq)(l + q - bq) - (1 + q - o)(l + g - b)q l 
Similarly B w and T w can be written as 



(5.7) 



and 



B = feg-/ 2 (l-g)(l + g) 

w (1 + g - ag)(l + g - 6g) - (1 + g - a)(l + g - %<" 1 ' ' 

(1 + 9) (? f (l-6 + g) + (6-l) g-1) 

^ = fa™ (a-l-g) (6-l-g)-((a-l)g-l) ((6-1) g-1))' (5 ' 9) 
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5.2. Limit w — > oo of the generating functions 
In the limit that w — > oo we have, for |g| < 1, 

L w {y/q/{l + g), a, 6) - ^ + q) , (5.10) 

1 + q — aq 

^M/(l + g ),a,6)^ (1 + ^ 1 ( t+ g -ag) (5 ' U) 

and B w (,/q/(l + q),a,b)->0. (5.12) 

We note that substituting z = z(q) back into the limits of L w and T w gives the generating 
functions of loops and tails adsorbing in a half-plane ( equations 13 . 61 and 13 . 7|) as required. 

5. 3. Singularities of the generating function of loops 

Since each of the models has the same denominator polynomial this confirms that all 
three models have the same limiting free energy. To determine the free energy we 
examine the dominant singularity of the denominator polynomial. 

The dominant singularity q c is the zero of the denominator of L w (q/(1 + q),a,b) 
that has the minimal value of \z\, and hence is the solution of 

„ = (l + q-aq)(l + q-bq) 
q (l + q-a)(l + q-b) V ; 

with this property. 

Before discussing the cases in which we can find the singularity analytically, we first 
consider some examples where we locate the singularity numerically. To begin let us 
consider the symmetric case a = b as this was also considered by DiMarzio and Rubin 
(1971). In figure El we plot K /oop (20; a, a) with the half-plane thermodynamic limit for 
comparison. For larger widths it is clear that the free energy converges to the half-plane 
curve as one might expect. Let us now consider weighting the vertices in the top wall 
such that they would strongly attract vertices of the polymer - one may expect this to 
hamper the convergence to the half-plane curve. Hence we consider 6 = 3 and again for 
width w = 20 plot the free energy K loop (20; a, 3): this can be found in figure El It can 
now be seen that for a < 3 the free energy is converging to something close to 0.75 and 
not the half-plane curve which is log(2) for a < 2. Physically, we can understand that a 
long polymer in a slit with a highly attractive top wall will stick to that wall while the 
attraction of the bottom wall is lesser in magnitude. However, this points to the fact 
that this large width limit for infinitely long polymers does not demonstrate the same 
physics as the half-plane case. 

6. Special parameter values of the generating function 



The solutions of the equation (|5.13|) cannot be written in closed form for all values of 
a and b. However, there are particular values for which this equation simplifies. When 
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Figure 5. A plot of the free energy for width w — 20 when a — b is shown as the 
diamond markers. The solid curve is the half-plane solution. 



a = 1 or 2 and 6 = 1 or 2, and when 

1 — a— equivalently ab = a + b (6.1) 

l — o 

the right-hand side of equation (J5.13j) simplifies. In these cases the locations of all the 
singularities of the generating functions can be expressed in closed form for all widths 
w. It is worth considering the generating function directly as some cancellations occur 
at these special points. 

For the special values mentioned the generating functions L w and B w have 
particularly simple forms, while the forms for T w are slightly more complicated. We 
give the explicit form for L w at each of the points below: 



when a = b = 1 



L - = — r^q~^ — ' ( } 



when a = 2 and 6=1 



Lw ~ (l-g)(l + r +l)' ^ 



when a = 1 and 6 = 2 



(l + g )(l + g ") 
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Figure 6. A plot of the free energy for width w = 20 when b — 3 is shown as the 
diamond markers. The solid curve is the half-plane solution. 



when a = b 



2 

L,. 



on the curve ab 
L,, 



{l + g){l + q w ) 

(l-g)(l-g™)' 

a + b 

(1 + q)((l + q-a) + (l + q- aq)q^ 



(6.5) 



(6.6) 



(1 + q - o)(l + q - aq)(l - q w ) 
The singularities in z of L w , B w and T w may be found from the singularities in q by 
the mapping ()5.1jl . The singularities of L w in q occur either for \q\ = 1 or for q G M + . 
The transformation (|5.1|) maps g G {|g| = 1} U 1R + to z 6 M + . We note that if q = e %e 
then 1/z — 2cos(#/2). We concentrate on the dominant singularity z c for which \z\ is 
minimal. 

Using the expressions for L w given above we find the zeros of the denominators 
(that do not cancel with a zero of the numerator) and choose the one that minimises 
\z\. This gives us q c and hence z c . From these we find the free energy k at the special 
points listed above: 

• when a = b = 1 the zero is q c = exp(2iri/(w + 2)) and hence 



k(w; 1, 1) 



log (2cos(tt/(w + 2))); 



(6.7) 
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• when a = 1, b = 2 and when a = 2, b = 1 the zero is q c = exp(iri/(w + 1)) and 
hence 

k(w; 2, 1) = k(w; 1, 2) = log (2 cos(vr/2(w + 1))) ; (6.8) 

• when a = b = 2 the zero is q c = 1 and hence 

«(w;2,2)=log(2); (6.9) 

• when ab = a + b and a > b the zero is q c = (a — 1) _1 while when ah = a + b and 
b > a the zero q c = {b — 1) _1 = (a — 1). Hence for a > 1 

a, a/(a - 1)) = log ( ^—^ J • ( 6 - 10 ) 

In a similar way, one may write down expressions for the locations of all the singularities 
of the generating functions. 

It is instructive to consider the free energy (see equation ([6.10)) ) along the special 
curve given by ab = a + b. We first note that it is independent of the width of the slit. 
It is also clearly different at all values of a < 2 to the half-plane free energy (which is 
log(2)). This re-iterates that the large width limit of the slit problem is not necessarily 
described by the half-plane results. 

In order to discuss the effective force between the walls for large w we need to 
calculate the large w asymptotics of the free energy. Since we have written down the 
closed form solution for all w of k(w; a, b) at some special points it is advantageous to 
consider the large w asymptotics of these first. As noticed above, on the curve ab = a + b 
the k is independent of w and so this includes the special point (a, b) = (2, 2). The large 
w asymptotics for k for the remaining special points are calculated by expanding in 
inverse powers of w. They are 

• when a = b = 1 

K(w,l,l)=log(2) 

• and when (a, b) = (1, 2) or (2, 1) 

k(w, 1,2) = k(w,2, 1) = log(2) 

7. Asymptotics for large widths at general parameter values 

As we have noticed the large w limit is dependent on the order of limits w —>■ oo and 
n — > oo. As such we can now consider the large w asymptotics of the free energy (which 
itself is first defined via ann->oo limit). As mentioned above, we are interested in the 
forces between the walls for large w which also involves large w asymptotics. However, 
we are unable to find a closed form solution for K,(w,a, b) for arbitrary a, b and w. We 
can, however, find asymptotic expressions for k for arbitrary a, b and large w without 
the need for such an expression. 

For general a and b the analysis depends on the region of parameter space and 
the plane naturally breaks up into various regions. These will become different phases 



vr 2 2tt 2 727t 2 + 7t 4 . s , , 

2^ + ^-^2^ + ^' ^ 



7T 2 7T 2 727T 2 + 7T 4 
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when considering the limit w — > oo. The analysis proceeds in a self-consistent manner 
perturbing around the solution obtained from the special points in descending powers 
of w. 

For a, 6 < 2 we can find the asymptotics of q c as a function of w, 

which, mapping back to the z variable, gives 

1 7T 2 7i 2 (ab-a-b) /5tt 4 37r 2 (a6 - a - 6) 2 \ 1 _ 5 
Zc ~ 2 + 4^ + (2 - a) (2 - + V 48~ + (2 - a) 2 (2 - 6) 2 ) ^ + ° (w } (7 ' 2) 
and so 

i / x 7T 2 2n 2 (ab — a — b) I V 4 6(a6 — a — 6) 2 \ 1 =-. .„ , 

« = log 2 - — = - - ^ rv-^ ~ ttt + TTT x,/o IvT " 5 )-(7-3) 



For 6 = 2 and a < 2 we have 

I + 

2 - a 



exp ( 7u/ ( w + ) + 0(w~ A ) ) (7.4) 



which gives 

1 7T 2 vr 2 a /5vr 4 3n 2 a 2 \ I 



and so 



2 + 16w 2 8(2-a)w 3 + ^768 + 16(2-a) 2 J +0(W 5) (T ' 5) 



7T 7T 2 a / 7T 4 37T 2 a 2 \ 1 r, . 

K = Og(2) + — r— j h T77 —7 + O(w^). (7.6) 

SV ; 8w 2 4(2 -a)w 3 \192 8(2 - a) 2 / w 4 K ' K ' 

When a = 2 and b < 2 we obtain a similar expression where b replaces a. 

When a or b is greater than 2 the asymptotic form changes. In particular the right- 
hand side of equation (|5.13|) is satisfied by either q — I /(a — 1) or q — 1/(6—1) and in 
the limit as w — > oo the left-hand side of the equation is zero. Hence we expand about 
the point q — l/(a — 1) and find that the next term in the expansion is exponential in 
w. More precisely, when a > 6 and a > 2 the asymptotic expansion of q c (w) as w — > oo 
is 

*-(^)- va:-v (^i)" + °H^r)- (? - 7) 

Transforming this back to the z variable gives 

a 2a(a-b)^fa^l \a - 1 J y \a-lj J K ' 

and so the free energy is 

For 6 > a and 6 > 2 one simply interchanges a and 6 in the above expressions to 
obtain the correct results. 
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Finally, when a = b > 2 the asymptotic expansion of q c (w) is 
and transforming this back to the z variable gives 



and 



a 2ay/a — 1 \a — l J \ \a — 1 



2 / i \ w/2 



«=n^) + fe$u^j +o in^rj j- (7i2) 

8. The infinitely wide slit 

We now come to consider the limit of large slits for infinite length polymers. Even 
though there isn't a closed form solution for finite widths of the equation (J5.13)) we can 
simply use the large w asymptotics above to deduce that 

log(2) if a, b < 2 

K mf ~ slu {a, b) = lim k(w; a,b) = { log (7=1) if a > 2 and a > 6 (g.l) 

log ( ^= J otherwise. 

One immediately sees that the free energy depends on both a and b rather than only 
on a which reflects the difference from the half-plane result f equation 13. 9j) noted earlier. 
In fact, the free energy is symmetric in a and b which reflects the observation that 
for infinitely long walks the end-points are irrelevant. For example, as in the case we 
considered earlier in figure |H1 with b > 2 the free energy for a < b is log (^7=) (which 
is approximately 0.752 for 6 = 3) while it is log( -^== ) for a > b. The half-plane phase 
diagram contains a desorbed phase for a < 2 and an adsorbed phase for a > 2. The 
transition between them is second order with a jump in the specific heat. Using (|8.1j) 
one can deduce for the infinite slit that there are 3 phases: one where the polymer is 
desorbed for a, b < 2; a phase where the polymer is adsorbed onto the bottom wall for 
a > 2 with a > b; and a phase where the polymer is adsorbed onto the top wall for 
b > 2 with b > a. This is illustrated in figure [3 The low temperature adsorbed phases 
are characterised by the order parameter of the thermodynamic density of visits to the 
bottom and top walls respectively. There are 3 phase transition lines. The first two are 
given by b = 2 for < a < 2 and a = 2 for < b < 2. These lines separate the desorbed 
phase from the adsorbed phases and are lines of second order transitions of the same 
nature as the one found in the half-plane model. There is also a first order transition 
for a = b > 2 where the density of visits to each of the walls jumps discontinuously on 
crossing the boundary non-tangentially. 
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b ♦ a = b 



b r =2 




a r = 2 



Figure 7. Phase diagram of the infinite strip. There are 3 phases: desorbed (des), 
adsorbed onto the bottom wall (ads bottom) and adsorbed onto the top (ads top). 



9. Forces between the walls. 

Let us define the effective force between the walls induced by the polymer as 

dn{w) 1 -dz c {w) 
= -r, = — 7— r — « • (9.1) 

aw z c (w) ow 

Before examining the general case for large w it is worth considering the general case for 
small w numerically and at the special points for all w where it can be found exactly. 
To begin, at the special points we can calculate the induced force exactly since we know 
the value of k(w) exactly. We note that since k(w) is monotonic at each of these points 
the induced force has one sign for all w. In particular, 

• when a = b = 1 then l/z c = 2 cos(7r/ (w + 2)) and so 

T = -^— 2 tan(vr/( W + 2)) = - 3 - **L + 0(^ 5 ) (9.2) 

which is positive and hence repulsive. 

• When (a, b) = (2, 1), (1, 2) then l/z c = 2 cos(vr/2(w + 1)) and so 

3tt 2 



Jr = 2(w 7T +iy tan{n/2{w + 1)) = 

which is positive and hence repulsive. 

When oh = a + b including a = b = 2 then l/z c 

F = 0. 



it 



4w 3 4w 4 



+ 0(w- 



(9.3) 



/a-l 



/b-1 



and so 



(9.4) 



If one considers more general values of a and b numerically one observes that in each 
region of the phase plane k(w) is monotonic and hence the force takes on a unique sign 
at each value of a and b. At general (a, b) we are able to find the asymptotic force for 
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F=7-z- 7T« + (9-6) 



large w by using the large w asymptotic expression for k(w; a, b) which we calculated 
above. 

Using the asymptotic expressions for k found in section in the (a, 6)-plane we 
obtain the asymptotics for the force. 

• For a, b < 2 

^_ 7r 2 6n 2 (ab — a — b) e. . . 

^" + (2- a) (2- 6 K +Q ''" » < 9 ' 5 » 

which is positive and hence repulsive. 

• For b = 2, a < 2 

7T 2 3vr 2 q 
4^3 _ 4(2 - a)w 4 
which is positive and hence repulsive. 
For a > b and a > 2 

(q-2) 2 (q&-q-6)log(q-l) / 1 \" // 1 \ 2w \ 
• F - 2(a-l)(a-6) U^I J + ° ^^TJ f ' U ' ] 

For b > a and b > 2 one simply interchanges a and b in the above expressions to 
obtain the required result. 

• At b = a > 2 

(q-2) 2 log(q- 1) / 1 \ w/2 // 1 \"\ 
J7= - 4(a-l) (^l) + °((^t) J' (9 ' 8) 
This is negative for all a > 2 and so the induced force is attractive. 

The regions of the plane which gave different asymptotic expressions for k and hence 
different phases for the infinite slit clearly also give different force behaviours. For the 
square < a, b < 2 the force is repulsive and decays as a power law (ie it is long-ranged) 
while outside this square the force decays exponentially and so is short-ranged. This 
change coincides with the phase boundary of the infinite slit phase diagram. However, 
the special curve ah = a + b is a line of zero force across which the force, while short- 
ranged on either side (except at (a, b) = (2, 2)), changes sign. Hence this curve separates 
regions where the force is attractive (to the right of the curve) and repulsive to the left 
of the curve. The line a = b for a > 2 is also special and, while the force is always short- 
ranged and attractive, the range of the force on the line is discontinuous and twice the 
size on this line than close by. All these features lead us to us a force diagram that 
encapsulates these features. This diagram is given in figure |H1 



10. Discussion 



We have solved and analysed in the limit of infinite length three types of polymer 
configuration, being loops, bridges and tails, in a two-dimensional slit geometry. The 
types of configurations considered are directed walks based on Dyck paths and the 
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short range 

repulsion 
' 

2 a 

Figure 8. A diagram of the regions of different types of effective force between the 
walls of a slit. Short range behaviour refers to exponential decay of the force with slit 
width while long range refers to a power law decay. The zero force curve is the special 
curve given by ab = a + b. On the dashed line there is a singular change of behaviour 
of the force. 

walks interact with both walls of the slit. We have found the exact generating functions 
at arbitrary width. We have calculated the free energy exactly at various points for 
arbitrary width, and, importantly, asymptotically (for large widths) for arbitrary a and 
b. This has allowed us to map out the phase diagram for infinite slits: for finite width the 
free energy is an analytic function of the Boltzmann weights a and b. This phase diagram 
is different to that obtained in the half-plane even though the generating functions for 
the half-plane are a formal limit of the finite width generating functions. This arises 
because the order of the infinite width and infinite polymer length limits (both of which 
are needed to see a phase transition) are not interchangeable. 

From the large width asymptotics we have mapped regions of the plane where 
the free energy decreases or increases with increasing width and, also, whether this 
happens with an algebraic decay or exponential decay. Using this we have delineated 
regions where the induced force between the walls is short-ranged or long-ranged and 
whether it is attractive or repulsive. Three types of behaviour occur: long-ranged 
repulsive, short-ranged repulsive and short-ranged attractive. Applying such a model 
to colloidal dispersions implies that the regions where the force is long-ranged and 
repulsive support steric stabilisation while the regions where the force is short-ranged 
and attractive promote sensitized flocculation. 

Given the curious interchange of limits phenomenon and from the point of view 
of the application of this theory to colloidal dispersions it would be interesting to 
investigate the finite polymer length cases both analytically and numerically with a 
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view of searching for a scaling theory for the crossover from half-plane to slit type 
behaviour. 

The model considered here is a directed walk in two dimensions. Although we 
expect the general features (such as the phase diagrams sketched in Figures 7 and 8) to 
be the same for self-avoiding walk models in both two and three dimensions, such models 
are beyond the reach of analytic treatments. It would be interesting to investigate the 
detailed behaviour of such models by Monte Carlo methods or by exact enumeration 
coupled with series analysis techniques. 
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